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摘 要 : 基于 印度 河流 域 及 周围 54 个 地 面 气象 站 气温 、 降 水 资料 ,结合 CRU 气温 和 GPCC 降水 全 球 
格 点 化 陆 面 再 分 析 资 料 ,通过 插值 构建 了 一 套 0.5° x 0. 5? 2-3 XE 1980 一 2016 年 逐 月 格 点 数据 集 。 
采用 Thornthwaite 方法 计算 了 洪 在 蒸 散 发 ,基于 标准 化 降水 蒸 散 指数 (SPE71) ,探讨 了 印度 河流 域 气 
候 变 化 及 于 旱 演 变 特征 。 结 果 表 明 :(1)1980 一 2016 年 ,印度 河流 域 年 平均 气温 以 0.30 C- 
(10a) ”的 速率 呈 显 著 上 升 趋势 ,21 世纪 初 增 温 幅度 最 大 ;和 干 季 (11 月 ~ 次 年 4 月 ) 升 温 速率 较 快 ， 
达 0.36 © - (10 a)”, 湿 季 (5 ~10 月 ) 增 速 0.25 C+ (10 a)”。 年 降水 量 呈 现 少雨 一 多雨 一 少 
雨 一 乡 雨 年 代 际 振荡 。 伴 随 着 持续 升温 ,年 和 各 季 的 潜在 蒸发 量 增 加 显著 。 干 季 干 旱 频率 较 多 ,但 
温 季 于 旱 强 度 高 ,各 季 干 旱 频率 与 降水 呈现 较 一 致 的 年 代 际 波动 ;干旱 的 影响 面积 在 干 季 呈 现 微弱 
地 增加 趋势 , 湿 季 却 略 有 减少 趋势 。(2) 空 间 上 , 除 西 北 局 部 ,流域 其 他 区 域 的 年 和 季 平 均 气 温 、 潜 
在 蒸发 量 增加 趋势 显著 , 均 达 到 95% 置信 水 平 。 其 中 南部 平原 和 东北 山区 升温 幅度 较 高 ,南部 平 
原 区 潜在 蒸发 量 增 加 也 较 大 。 新 德里 到 喀布尔 的 东南 至 西北 带 状 区 域 的 年 和 湿 季 降水 量 ,以 及 喀 


南 局 部 山区 的 湿 季 呈现 显著 的 干 
增多 带 来 的 不 利 影响 。 
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2006—2015 年 ,全 球 地 表 平 均 温度 较 1850 一 
1900 年 升 高 了 0.87 % (0.75 ~ 0.99 ©) ,2017 年 
升 高 1 "TI 。 全 球 变 暖 背景 下 ,极端 天 气 气 候 事件 
越 来 越 频 繁 ,干旱 强度 .频率 和 持续 时 间 不 断 增加 ， 
对 水 资源 .生态 系统 和 社会 经 济 产 生 严 重 的 负面 影 
mg’? *) 。 气 候 变化 及 干旱 演变 特征 的 研究 ,可 以 为 
灾害 风险 管理 减轻 干旱 对 社会 经 济 可 能 造成 的 影 
响 和 制定 适应 干旱 灾害 的 战略 提供 科学 依据 。 

南亚 印度 河流 域 是 世界 上 面积 最 大 的 流域 之 
一 ,中 下 游 遍 布 着 世界 最 大 的 灌溉 网 络 ,是 农业 经 济 
为 主导 的 国际 河流 ,大 约 2.15 x10* 人 口 的 生计 直 
接 或 间接 地 依赖 于 农业 和 其 他 与 流域 水 资源 相关 职 
Jp? 75 。 流 域 大 部 属于 干旱 、 半 干旱 气候 ,降水 稀少 
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布尔 周围 地 区 的 干 季 降 水 量 呈 显著 增加 趋势 。 东 南平 原 区 和 东北 局 部 山区 的 干 季 , 以 及 东北 和 西 
旱 化 态势 ,需要 加 强 防 灾 减 灾 的 意识 并 采取 相应 措施 ,以 规避 干旱 


气温 ; 降水 ; ABA; 干旱 ; 时 空 变化 ; 印度 河流 域 


且 变 率 大 ,生态 环境 脆弱 ,对 全 球 气 候 变 化 十 分 敏 
感 "”-” 。 频 繁 发 生 的 干旱 ,对 区 域 农 业 生产 造成 了 
严重 的 影响 ,其 至 威胁 到 粮食 安全 ""“ 。 另 外 ,由 
于 气候 变化 带 来 的 水 资源 不 稳定 性 ,使 位 于 干旱 半 
干旱 区 的 这 类 国际 性 河流 面临 着 跨 境 水 资源 安全 及 
纠纷 问题 -5 。 已 有 研究 对 印度 河流 域 分 区 的 气 
温 .降水 .干旱 等 的 演变 进行 了 分 析 , 发 现 1967 一 
2005 年 印度 河流 域 上 游 以 及 1971 一 2010 年 中 、 下 
游 气温 呈 上 升 趋势 ,降水 却 没有 任何 显著 的 变化 趋 
416-17) ;1960 一 2013 年 和 1980—2016 年 印度 河流 
域 巴基斯坦 境内 变 暖 趋势 明显 , 尤 以 南部 地 区 升温 
速率 最 快 开 -8 ;1980 一 2016 年 巴基斯坦 北部 和 中 
部 降水 量 呈 增加 趋势 , 且 在 北部 高 海拔 山区 变 率 较 


基金 项 目 : 科技 部 基础 资源 调查 专项 (2018FY100501 ) ;国家 自然 科学 基金 委 和 巴基斯坦 科学 基金 会 合作 项 目 (41661144027) ;国家 自然 


科学 基金 (41671211 ) 资助 


作者 简介 : 赵 建 婷 (1993 - ) , 女 ,硕士 研究 生 , 研 究 方向 为 气候 变化 影响 评估 . E-mail :zjmuist_2017@ 163. com 
通讯 作者 : 姜 彤 , 男 , 研 究 员 ,研究 方向 为 气候 变化 风险 管理 . E-mail :jiangtong@ cma. gov. cn 


202006.00024v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


$4 d wy 


大 ,南面 和 西南 部 分 平原 区 却 呈 减少 趋势 ”|。 
1961—2007 年 印度 河流 域 西部 和 南部 频繁 遭受 中 
度 及 以 上 的 干旱 ,尤其 在 1998 一 2002 年 经 历 了 历史 
上 最 严重 的 干旱 ,导致 巴基斯坦 南 及 西南 省 份 遭 受 
严重 饥荒"”-” ;1986 一 2005 年 印度 河流 域 每 年 干 
旱 影响 面积 约 260. 40 x 10° km ,影响 人 口 4.32 x 10" 
人 3。 此 外 ,1982 4Æ 1987 年 、1996 年 .2000 年 、 
2002 年 和 2009 年 印度 极端 干旱 累积 影响 人 口 超过 
7.51 x105i] 。 上 述 研究 多 采用 巴基斯坦 境内 测 站 
资料 ,以 印度 河流 域 局 部 区 域 为 研究 单元 ,未 关注 流 
域内 其 他 区 域 的 气候 变化 特征 。 对 干旱 研究 主要 基 
于 降水 ,未 考虑 气温 .蒸发 等 要 素 的 影响 , 且 研 究 时 
段 以 年 尺度 为 主 , 对 季节 性 干旱 的 研究 尚 有 不 足 。 
本 研究 收集 整理 了 印度 河流 域 及 周围 地 区 地 面 
观测 资料 ,结合 CRU 气温 和 GPCC 降水 再 分 析 资 
料 ,构建 了 覆盖 整个 印度 河流 域 的 0.5° x0.5° 分 辨 
率 数据 集 。 采 用 标准 化 降水 落 散 指数 (SPE1)， 
探讨 了 印度 河流 域 气候 变化 和 干旱 演变 特征 ,以 期 
为 印度 河流 域 应 对 气候 变化 和 抗旱 减灾 等 措施 的 制 
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定 提供 科学 依据 ,为 国际 河流 水 资源 管理 提供 重要 
借鉴 。 


1 研究 区 概况 和 方法 


1.1 研究 区 概况 

印度 河流 域 位 于 南亚 , 介 于 22° ~ 38°N 和 
65° ~ 87"E 之 间 ,面积 约 为 110 x 10* km ,由 阿富汗 
(796) , P E (1196) , Ep BE (27% ) 和 巴基斯坦 
(55% )4 个 国家 共同 拥有 "“…] 。 干 流 印度 河 (全 长 
3 180 km) 发 源 于 中 国 青藏 高 原 , 流 经 巴基斯坦 , 注 
入 阿拉 伯 海 ;喀布尔 河 起 源 于 阿富汗 , 流 经 喀布尔 市 
进入 巴基斯坦 境内 并 汇 入 印度 河 , 是 印度 河 右 侧 的 
主要 文 流 ;而 杰 赫 和 勒 姆 河 . 奇 纳 布 河 . 拉 维 河和 萨 网 
鲁 伊 河 为 向 西 流入 流域 的 主要 河流 ”9 。 流 域 地 形 
起 伏 差异 明显 ,北部 以 山地 为 主 ,多 数 地 区 海拔 在 
7 000 m 以 上 ,最 高 点 达到 8 569 m; 南 部 和 阿拉 伯 
海 沿岸 地 区 为 平坦 的 农业 平原 ,海拔 最 低 点 为 
-44 m^?! (图 1)。 
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Fig.1 Location of hydrological and meteorological stations in the Indus River Basin 


1.2 资料 
1.2.1 气象 观测 资料 ”收集 整理 了 印度 河流 域 及 


周边 区 域 167 个 地 面 气象 站 观测 资料 ,包括 巴 基 斯 
坦 气 象 局 提供 的 50 个 地 面 观 测 站 逐日 气温 和 降水 ， 
其 中 45 个 站 点 资料 时 段 为 1961 一 2013 年 ,5 个 站 


点 资料 起 始 时 间 略 偏 晚 (分 别 为 1963 年 ,1968 年 ， 
1971 年 , 1974 年 , 1975 年 ); 110 个 国际 交换 站 
1951 一 2014 年 逐日 气温 ,其 中 80 个 位 于 印度 ,23 个 
位 于 阿富汗 ,7 个 位 于 塔吉克 斯 坦 ;7 个 中 国 气象 局 
1951 一 2017 年 逐日 气温 和 降水 。 


赵 建 婷 等 :印度 河流 域 气温 降水 ,蒸发 及 干旱 变化 特征 


考虑 到 记录 的 完整 性 和 可 比 性 , 当 测 站 缺失 年 
份 超过 5% ,该 站 未 参加 计算 ; 当年 内 缺失 超过 5% , 
该 年 资料 未 参加 计算 。 经 过 严格 的 气候 界限 值 
检查 和 气候 极 值 检查 1 ,研究 保留 了 38 个 巴 基 斯 
坦 气象 站 ,1 个 塔吉克 斯 坦 和 8 个 印度 国际 交换 站 ， 
以 及 7 个 中 国 气象 观测 站 , 共 54 个 测 站 的 观测 资料 
(图 1)。 由 于 53% 的 观测 站 资料 在 1960s 和 1970s 
缺失 天 数 占 总 缺失 天 数 的 比例 超过 80% , 而 1980 
年 后 记录 比较 完整 ,文中 分 析 了 1980 年 以 来 的 印度 
河流 域 气 候 变化 特征 。 流 域 降水 受到 季风 系统 的 控 
制 , 具 有 明显 的 干 湿 季 节 ,5 ~ 10 月 为 湿 季 ,11 月 ~ 
次 年 4 月 为 干 季 。 
1.2.2 水 文 观 测 资 料 径流 数据 来 源 于 巴基斯坦 
水 电 开 发 局 (Water and Power Development Authority , 
WAPDA) ,包括 印度 河流 域 上 游 Nowshehra , Tarbela 
和 Mangla 3 个 水 文 站 1980 一 2009 年 逐日 流量 。Now- 
shehra 水 文 站 是 喀布尔 河水 文 站 ,控制 流域 面积 约 
88 676 km" ,流域 内 冰川 覆盖 面积 约 2% 77 ; 塔 贝 拉 
流域 Tarbela 水 文 站 控制 流域 面积 约 2 x 10° km? ,其 
中 约 12% 的 流域 面积 被 冰川 覆盖 ; Mangla 是 东部 
文 流 杰 赫 勒 姆 河水 文 站 ,流域 面积 约 为 34 037 
km) 。 印 度 河流 域 水 文 站 空间 分 布 见 图 1 所 示 。 
1.2.3 再 分 析 资料 英国 东 英 吉利 大 学 气候 研究 
中 心 (Climatic Research Unit , CRU ) 的 气温 和 美国 全 
球 降水 气候 中 心 (Global Precipitation Climatology 
Centre ,GPCC ) 的 降水 资料 是 当前 应 用 最 广泛 的 两 
套 全 球 格 点 化 陆 面 资 料 ” 2016 年 发 布 的 CRU 
TS v4.01(CRU Time Series Version 4.01) 气 温 资 料 , 
数据 长 度 为 1901 一 2016 年 ,空间 分 辨 率 为 0.5? x 
0.5°;2018 年 发 布 的 GPCC FDP v8 (Full Data Rean- 
alysis Version 8) 降水 资料 ,时 间 尺 度 为 1891—2016 
年 ,空间 分 辨 率 为 0.5" x 0.5"。 本 文选 取 1980 一 
2016 年 的 CRU 逐 月 气温 和 GPCC 逐 月 降水 以 辅助 
观测 资料 的 格 网 化 。 
1.3 方法 
1.3.1 潜在 基 散 发 计算 ”潜在 蒸 散发 是 指 水 分 供 
给 不 受 限时 ,下 垫 面 可 能 达到 的 最 大 闵 发 量 ” 。 文 
中 选择 计算 方案 最 为 方便 、 数 据 需 求 最 少 的 Thorn- 
thwaite 方法 ,分 析 了 潜在 蒸 散 发 演变 及 干旱 时 空 
分 布 。 其 计算 公式 如 下 : 

0 T,<1 °C 


107, 
PET, = 1.6d[—+]* x10 


1C <T,<26.5C (1) 


a,*a,T,-a,T^ T,>26.5 C 


式 中 :d 为 每 月 的 天 数 除 以 30;a, = - 415. 854 7; 
a, 232.244 1;a, =0. 432 5, T, 为 月 平均 温度 ;a = 
0. 492 39 + 1. 792(10) 71 - 7. 71 (10) ^P +6.75 
(10) "P f= Xi XA AREG i = [7/5] 98 
月 平均 加 热 指数 ;7 为 第 i 个 月 的 月 平均 温度 。 
1.3.2 干旱 指教 ”和 常用 的 干旱 指数 有 帕 尔 默 干旱 
指数 (Palmer Drought Severity Index , PDS/) D5] 标准 
4k WE Jk 48 BL ( Standardized Precipitation Index, 
SPI) °°) 标准 化 降水 蒸 散 指数 ( Standardized Precipi- 
tation Evapotranspiration Index , SPET) [24,37] 和 标准 化 
径流 指数 (Standardized Discharge Index , SDI) [38] A 
其 中 ,SPEI 综合 考虑 气温 、 降 水 和 莹 散发 对 干旱 的 
影响 ,适用 于 不 同时 间 尺 度 旱 情 监测 和 评估 7^ 
SDI 利用 观测 流量 ,用 工分 布 描述 径流 量 的 变化 并 

经 正 态 标 准 化 求 得 ,能 客观 反映 径流 干枯 变化 。 
SPEI 的 计算 ,首先 分 析 逐 月 降水 量 与 蒸 散 量 的 
PAE ,并 采用 三 参数 的 log-logistic 概率 分 布 对 差 值 

序列 进行 正 态 化 ,再 对 累积 概率 密度 进行 标准 化 : 
D, = P, - PET, (2) 


F(x) = [1+ (2—)6] (3) 


2 
Co 十 C1W 十 C220 


l+dw+d,w, +d,w, 


SPEI =w - 


(w= /-2In(p) , px0.5) (4) 


2 
Co CW t C310 


1+d,wt+d,w, + d4w, E 


(w= /-2In(1 =p), p>0.5) (5) 
式 中 :PP; 为 降水 量 ,PE7 AN ZR it x 表示 自 变 量 D, 
值 ,a 为 尺度 参数 ,B 为 形状 参数 ,y 为 位 置 参 数 ;p 
为 累积 概率 ,常数 项 cv =2. 515 517,c, =0. 802 853, 
c, 20. 010 328 ,d, =1. 432 788 ,d =0. 189 269 ‚d; = 
0.001 308'*! 。 
SDI 计算 与 SPEARS]; 


SPEI = 


E t- (et - c )t +c, 
Beh? aed oa edad (6) 


In u 
G(x)’ 


(7) 


CE NL x, (x 
GG) | #dx,(x>0) (8) 
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式 中 :5 为 概率 密度 正 负 系数 ;x 为 径流 样本 值 ;C 
(x%) 为 与 T 函数 相关 的 径流 分 布 概率 ; 当 G(x) > 
0.5 时 为 正 ,G(x) «0. 5 时 为 负 。 

本 研究 选择 6 个 月 尺度 SPEI (SPE1_6) 表 征 干 


旱 , 并 采用 6 个 月 尺度 SDI CSDI 6) 进行 比较 分 析 。 
以 SPEI < -1 的 格 点 占 总 格 点 数 的 百分比 来 评价 干 
旱 影响 面积 的 大 小 及 其 严重 程度 。SPA1 和 SDI 的 
干旱 等 级 划分 见 表 1 所 示 A, 


表 1 SPEI 和 SDI 指数 干旱 等 级 划分 
Tab.1 SPEI- and SDI- based drought classifications 


SPEI/SDI < -2.0 [-2.0,-1.5) [-1.5,-1.0) [-1.0,1.0) [1.0,1.5) [1.5,2.0) 22.0 
FHER 极端 干旱 重度 干旱 中 度 干 旱 正常 PRE EE HKI 


1.3.3 数据 格 网 化 方法 ”以 CRU 多 年 平均 月 气温 
和 GPCC 多 年 平均 月 降水 量 作为 月 气温 、 降 水 的 平 
均 气 候 态 ,采用 反 距离 加 权 方 法 (IDW ) 将 平均 气候 
态 插 值 到 地 面 观测 站 ,获取 与 地 面 观测 值 的 距 平 。 
其 中 ,气温 距 平 为 气候 场 与 观测 值 的 差 值 , 降 水 距 平 
为 气候 场 与 观测 值 的 比值 。 将 此 距 平 采用 IDW 
方法 ,插值 并 又 加 到 气候 态 格 点 (气温 相 加 ,降水 相 
HE) ,获得 0.5° x0.5° 分 辨 率 的 格 网 数据 集 ， 以 拟 补 
地 面 观测 数据 的 空间 不 均匀 性 。 


2 结果 与 分 析 


2.1 气温 时 空 演变 

1980—2016 年 印度 河流 域 年 平均 气温 以 0. 3 
C+ (10 a) ! 的 速率 量 显著 上 升 趋势 ,达到 99% fri 
度 水 平 。 其 中 ,2016 年 达到 最 高 为 19.5 C ,1983 年 


20.0 [ 
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最 低 为 17.5 (图 2a)。 从 气温 距 平 的 9 a 滑动 平 
均 看 ,1996 年 后 气温 距 平 由 负 转 为 正 ,21 世纪 初 升 
温 幅 度 急剧 增 大 ,2000 年 后 ,气温 距 平 基 本 上 都 在 
0.24 % 以 上 ,2003 年 达到 最 大 距 平 0.43 % (图 
2b) 。 空 间 上 , 除 西北 局 部 山区 ,其 他 86. 6% 的 面积 
升温 趋势 显著 ,通过 95% 信和 度 检 验 。 升 温 幅 度 在 流 
域 南部 平原 和 东北 山区 较 高 ,为 0.48 ~ 0.59 "C - 
(10 a) ,西部 相对 较 低 ,0.04 ~ 0.15 CT: 
(10 a) (图 5a)。 

干 湿 季 平均 气温 上 升 趋势 同样 显著 ,但 增 速 略 
有 差异 。 干 季 升 温 速率 达 0.36 "C - (10 a) ,气温 
PEF 1998 年 后 由 负 转 为 正 ,2000s 中 期 气温 增加 最 
为 剧烈 , 自 2000s 后 期 以 来 升温 速率 略 有 减缓 ; 湿 季 
增 速 0.25 C+ (10 a) 5,1997 年 后 距 平 由 负 转 为 
IE , 2000s 前 期 升温 幅度 急剧 增 大 ,2000s 中 期 到 后 
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图 2 1980—2016 年 印度 河流 域 年 、 季 平均 气温 及 其 距 平 的 9 a 滑动 平均 
Fig.2 Variation of annual and seasonal temperature and 9 year moving average of temperature anomalies in 


the Indus River Basin during 1980 -2016 


赵 建 婷 等 :印度 河流 域 气温 .降水 ,蒸发 及 干旱 变化 特征 


期 略 有 组 减 , 自 2010 年 以 来 增幅 又 有 所 增加 (图 2c 
~2d)。 流 域 大 部 地 区 干 湿 季 气温 上 升 趋势 显著 。 
其 中 , 干 季 升 温 幅 度 在 东北 山区 及 东南 平原 区 相对 
较 高 ,西部 山区 相对 较 低 ; 湿 季 东北 山区 最 大 ,西北 
山区 西部 区 域 最 小 (图 略 )。 

2.2 降水 时 空 演变 

1980—2016 年 ,印度 河流 域 年 降水 量 呈 少雨 一 
多 雨 一 少雨 一 多 雨 年 代 际 振荡 ,1995 年 达到 最 大 值 
为 784.2 mm,2002 年 最 小 值 为 271.9 mm( 图 3a ~ 
3b)。 从 空间 变化 看 , 除 流域 北部 、 东 北 局 部 山区 ， 
其 他 区 域 降水 量 均 呈 增 加 趋势 ,在 新 德里 延伸 至 喀 
布尔 的 东南 西北 向 带 状 区 域 ( 占 流域 32. 6% ) 增加 
趋势 为 显著 (95% 信和 度 水 平 ) ,增幅 超过 60 mm - 
(10 a) (图 5$b) 。 

就 季节 而 言 , 干 季 1980s 中 期 和 1990s 末期 一 
2000s 中 期 表现 为 少雨 ,1980s 后 期 一 1990s 中 期 和 
自 2000s 后 期 以 来 呈 多 雨 ,多雨 期 和 少雨 期 最 大 距 
平 超过 20 mm; 湿 季 在 1980s 中 后 期 和 2000s 表现 
为 少雨 ,1990s 和 2010—2016 年 表现 为 多 雨 ,多 雨 期 
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显著 增加 趋势 (99% f E), XE 17. 9 mm - 
(10 a) ^! ,2010 年 达到 最 高 1 148.5 mm ,1983 年 最 
低 为 1 039.9 mm( 图 4a)。 从 潜在 蒸发 距 平 的 9 a 
滑动 平均 看 ,1997 年 后 潜在 莱 发 量 的 距 平 由 负 转 为 
正 ,21 世纪 初 增 速 急 剧 加 大 (图 4b)。 空 间 上 ,与 气 
温 变 化 趋势 相似 , 除 西 北 局 部 山区 ,流域 85.7 % 的 
地 区 增加 趋势 显著 (95% 信和 度 ) , 增 速 最 大 的 地 区 为 
东南 平原 区 域 ,为 31.5 ~39.0 mm * (10 a) ;西部 
山区 增幅 相对 较 小 ,为 0.5 ~8.2 mm - (10 a) (图 
36)s 

干 湿 季 潜 在 蒸发 量 分 别 以 9.1 mm + (10a) * 
和 8.9 mm * (10 a) 的 速率 呈 增 加 趋势 ,达到 99% 
的 信和 度 水 平 。 干 季 淤 在 蒸发 量 距 平 1998 年 后 由 负 
转 为 正 ,2000s 中 期 增加 最 为 剧烈 ,2000s 后 期 以 来 
增 速 略 有 减缓 ; 湿 季 1997 年 后 距 平 由 负 转 为 正 ， 
2000s 前 期 升幅 最 大 ,2000s 中 期 到 后 期 咯 有 组 减 ， 
自 2010 年 以 来 增幅 有 所 增加 (图 4e ~4d) FER 
在 蒸发 增加 趋势 在 流域 大 部 区 域 均 显 著 (95% 信 
JE) ,增加 幅度 干 季 在 中 部 和 东南 平原 区 , 湿 季 在 中 


和 人 少雨 期 最 大 距 平 超过 40 mm (图 3c ~3d) 。 空 间 
上 ,西北 喀布尔 局 部 山区 干 季 ,东南 新 德里 延伸 至 西 
北 喀布尔 的 带 状 区 域 湿 季 降 水 量 增加 显著 (95% 信 
BE) ; 北 、 东 北 和 西部 局 部 山区 干 湿 季 降水 量 ,以 及 
中 部 干 季 降 水 量 则 有 不 明显 的 减少 趋势 (图 略 ) o 
2.3 ”潜在 蒸发 量 时 空 演变 

伴随 着 显著 变 暖 ,印度 河流 域 年 潜在 蒸 散发 呈 
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部 、 东 部 平原 及 东北 局 部 山区 最 大 (图 略 )。 

2.4 干旱 时 空 演变 

2.4.1 SPEI 5 SDI 比较 为 验证 SPEI 对 印度 河 
流域 干旱 程度 的 描述 能 力 ,基于 流域 上 游 的 喀布尔 
子 流域 塔 贝 拉 子 流域 和 杰 赫 勒 姆 子 流域 的 3 个 水 
文 站 1980—2009 年 逐日 流量 计算 获得 SD1_6, 并 与 
相应 子 流 域 的 SP&k1_6 进 行 了 比较 ,由 图 6 可 见 ,2 
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图 3 1980—2016 年 印度 河流 域 年 季 降 水 量 及 其 距 平 的 9 a 滑动 平均 


Fig.3 Variation of annual and seasonal precipitation and 9-year moving average of precipitation anomalies in 


the Indus River Basin during 1980 -2016 
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图 4 1980—2016 年 印度 河流 域 年 . 季 潜 在 蒸发 量 及 其 距 平 的 9 a 滑动 平均 


Fig.4 Variation of annual and seasonal potential evapotranspiration and 9-year moving average of potential evapotranspiration 


anomalies in the Indus River Basin during 1980 — 2016 
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种 指数 能 较 一 致 地 反映 喀布尔 流域 . 塔 贝 拉 流域 和 
杰 赫 勒 姆 流域 干 湿 季 干 时 演变 特征 。 其 中 ,喀布尔 
流域 干 季 、 塔 贝 拉 流域 干 和 于 和 杰 赫 勒 姆 流域 干 湿 季 
的 相关 性 较 好 ,两 种 指数 的 相关 系数 达 0. 82.0. 58、 

0.87 和 0.42(99% 信 度 ) ,喀布尔 流域 湿 季 相关 系 
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TE: 白 点 表示 变化 趋势 通过 95% 信 度 检验 
图 5 1980 一 2016 年 印度 河流 域 年 气温 、 降 水 量 和 
潜在 蒸发 量变 化 趋势 的 空间 分 布 
Fig.5 Spatial distribution of trends of annual temperature, 


precipitation and potential evapotranspiration over 


the Indus River Basin during 1980-2016 


BOA 0. 38 (95% fii RE 
略 低 。 

从 SPEL6 逐年 变化 情况 看 , 干 季 ,1999 一 2001 
年 3 个 子 流域 都 出 现 连 年 中 度 及 以 上 干旱 ,1984 和 
1985 年 塔 贝 拉 流 域 和 杰 赫 勒 姆 流域 也 有 中 度 干 早 


) , 塔 贝 拉 流域 湿 季 相关 人 性 
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发 生 。 湿 季 ,1999 一 2001 年 喀布尔 流域 发 生 连 年 中 
度 及 以 上 干旱 .2001 年 塔 贝 拉 流 域 发 生 极端 干旱 、 
2000 年 杰 赫 勒 姆 流域 发 生 中 度 干 旱 , 与 SDI 指数 反 
映 的 干旱 年 份 和 等 级 有 着 较 好 的 对 应 状态 。 仅 个 别 
年 份 的 湿 季 ,SPEI1_6 反映 比 SDI_6 严重 的 干旱 程度 
(图 6)。 

2.4.2 干旱 时 宝 变化 ， 根 据 SPEI 6,1980—2016 
年 印度 河流 域 干 季 干旱 频率 较 高 ,37 a 里 共有 9 a 
发 生 了 中 度 及 以 上 的 干旱 ;但 从 干旱 等 级 看 , 湿 季 更 
容易 发 生 较 强 干旱 ,共有 4 次 重度 干旱 ,1 次 极端 干 
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图 6 1981—2009 年 印度 河上 游 喀 布尔 流域 , 塔 贝 拉 流域 和 杰 赫 勒 姆 流域 干 季 
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早 。 干 季 , 干 旱 主要 集中 在 1980s 中 期 和 1990s 后 
期 至 2000s。 其 中 ,后 段 干 旱 频率 多 且 等 级 较 高 ， 
1999—2008 年 间 , 共 有 7 a 达到 中 度 及 以 上 等 级 干 
旱 ,2000 年 .2001 年 和 2004 年 达到 重度 干旱 (图 
Ta). WZ ,1980s 一 2000s ,干旱 频率 不 断 增加 ,等 级 
不 断 加 重 。1990s 后 期 一 2000s 发 生 史 上 最 严重 的 
连年 干旱 ,1999 年 .2001 年 达到 重度 干旱 ,2002 年 
达到 极端 干旱 ;另外 ,2009 年 也 发 生 了 极端 干旱 (图 
Tb) 。 

就 干旱 面积 而 言 , FZ, 1980s 中 期 和 1990s 后 
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湿 季 SPEI 6 和 SDI_6 变化 


Fig.6 Comparison of dry and wet seasonal SPEI 6 and SDI 6 for the Kabul River UN Tarbela River Basin and Jhelum 
River Basin in the Upper Indus River Basin during 1981 — 2009 
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图 7 1980—2016 年 印度 河流 域 干 季 MEE SPEI 指数 时 间 变 化 


Fig.7 Variation of dry and wet seasonal SPEI in the Indus River Basin during 1980 - 2016 
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图 8 1980—2016 年 印度 河流 域 干 季 \ 湿 季 干 旱 面积 变化 


Fig.8 Variation of percentage of drought stricken areas over the Indus River Basin for wet and dry seasons during 1980 — 2016 
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图 9 1980—2016 年 印度 河流 域 干 季 、 湿 季 SPEI 变化 趋势 空间 分 布 
Fig.9 Spatial distribution trends of SPEI in the Indus River Basin for dry and wet seasons during 1980 —2016 


期 一 2000s 影响 面积 较 大 ; 1980s 前 期 和 1980s 后 
期 一 1990s 中 期 相对 较 小 ; 湿 季 ,1980s 及 2000s 以 
来 频繁 发 生 影响 面积 较 大 的 干旱 ,而 小 面积 的 干旱 
多 集中 发 生 于 1990s (图 8). 干 季 , 中 度 和 重度 干 
旱 均 呈 上 升 趋势 ,增长 率 分别 达 1.56 % - (10a) ^ 
和 0.12 % - (10 a) ^; WEE, 2 种 等 级 干旱 影响 面 
积 都 有 减少 趋势 ,减少 率 分 别 为 -0.31 96 - (10 
a) 和 -0.15 % - (10a) ^', 

干旱 强度 , 干 季 除 西北 局 部 山区 有 微弱 的 减弱 
趋势 外 ,流域 73.8% 的 面积 上 表现 出 不 同 程度 的 加 
剧 趋 势 。 其 中 ,流域 东南 平原 区 及 东北 局 部 山区 约 
20. 4% 的 面积 增加 趋势 显著 (95% 信和 度 ) (图 9a) ; 相 
反 , 湿 季 流 域 大 部 区 域 干旱 强度 则 呈 减 弱 态 势 , 仅 东 
北 和 西南 小 范围 山区 呈 加 剧 态势 (图 9b)。 


3 结论 


全 球 气候 变化 是 当今 各 国政 府 所 关注 的 焦点 问 
题 。 干 旱 作 为 最 频繁 发 生 的 自然 灾害 之 一 ,对 社会 、 


经 济 和 环境 产生 的 影响 重大 。 随 着 全 球 气温 持续 天 
高 ,印度 河流 域 气候 及 干旱 特征 也 发 生 相 应 的 变化 。 
本 文 基于 54 个 印度 河流 域 及 周边 区 域 地 面 气 象 站 
的 气温 和 降水 观测 资料 ,结合 CRU 气温 和 CPCC 降 
水 全 球 格 点 化 陆 面 再 分 析 资 料 ,构建 了 一 套 0.5° x 
0.5* 分 辩 率 逐 月 格 点 数据 集 。 采 用 Thomthwaite 方 
法 计算 了 潜在 蒸 散 发 量 ,基于 标准 化 降水 燕 散 指数 
(SPET) ,对 印度 河流 域 气候 变化 及 干旱 演变 特征 展 
开 研究 。 结 论 如 下 : 

(1) 1980 一 2016 年 ,印度 河流 域 年 平均 气温 以 
0.30 © + (10 a) ”的 速率 呈 显 著 上 升 趋势 ,21 世纪 
初 升温 幅度 急剧 增 大 ; 干 季 (11 月 ~ 次 年 4 月 ) 升 温 
速率 较 快 , 达 0.36 "C (10 a) ^ SEE (S ~10 H) 
增 速 0.25 "C + (10 a) -，。 除 西北 局 部 区 域 ,其 他 区 
域 的 年 和 季 气 温 上 升 趋 势 显著 (95% 的 信 度 水 平 )， 
其 中 ,南部 平原 和 东北 山区 升温 幅度 较 高 ,西北 局 部 
山区 相对 较 小 。 年 降水 呈 少 雨 一 多 雨 一 少雨 一 多 十 
年 代 际 振荡 。 东 南 新 德里 至 西北 喀布尔 流域 的 带 状 
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区 域 的 年 和 湿 季 降水 量 , 以 及 喀布尔 周围 区 域 干 季 
降水 量 呈 显著 增加 趋势 ,其 他 区 域 没有 明显 的 变化 
特征 。 伴 随 持续 升温 , 洪 在 蒸发 量 呈 明显 增加 趋势 ; 
空间 格局 与 气温 一 致 , 变 率 在 南部 平原 较 大 ,西北 局 
部 山区 较 小 。 

(2) 印度 河流 域 干 季 干 旱 频率 较 高 ,但 湿 季 干 
AGRE K, FÆ 1980s 中 期 和 1990s 后 期 一 2000s 
少雨 期 干旱 频率 最 高 , 且 干 时 影响 面 积 较 大 。 湿 季 
自 1980s 一 2000s ,干旱 频率 不 断 增 加 , 且 强 度 不 断 加 
重 ,1990s 后 期 一 2000s 少雨 时 期 ,发 生 连年 重度 及 
以 上 干旱 ,干旱 影响 面积 最 大 ; 干 季 干旱 影响 面积 呈 
现 出 增加 趋势 , 湿 季 则 有 减少 趋势 。 

(3) 气候 变 暖 TALE ZS Ae at b SS IY Be P LED 
度 河流 域 干 湿 季 干旱 强度 变化 的 空间 差异 较 大 。 F 
季 , 除 西北 局 部 山区 干旱 强度 有 减弱 趋势 外 ,其 他 区 
域 呈现 不 同 程度 的 加 剧 趋势 ,以 东南 平原 区 及 东北 
局 部 山区 趋势 显著 ; 湿 季 大 部 区 域 干旱 趋 于 减弱 ,但 
东北 和 西南 局 部 山区 呈 加 剧 态势 。 


4 讨论 


本 研究 基于 印度 河流 域 及 周围 地 面 气象 站 资料 
和 全 球 再 分 析 数 据 ,构建 一 套 1980 一 2016 年 0.5° x 
0.5° 分 辨 率 的 格 点 数据 集 。 根 据 本 研究 所 构建 的 气 
温 序列 ,印度 河流 域 年 平均 气温 在 南部 平原 区 超过 
25 C ,北部 高 山区 在 5 % 以 下 ,与 CRU 气温 产品 的 
差别 基本 控制 在 5 % ,两 者 空间 相关 系数 达到 
0. 99 ;而 在 高 海拔 和 复杂 地 形 过 渡 区 域 , 较 CRU 数 
据 偏 高 3 以上。 根据 本 研究 所 构建 的 降水 序列 ， 
印度 河流 域 年 降水 在 北部 山区 超过 1 300 mm, 而 南 
部 平原 .沙漠 区 却 不 到 200 mm ,与 GPCC 产品 在 流 
域 大 部 区 域 差 异 基本 控制 在 5% ,两 者 空间 相关 系 
数 为 0.98 ; 而 在 东北 部 山区 和 东部 局 部 地 区 差异 明 
显 ,其 中 东部 局 部 地 区 较 GPCC 偏 大 ,东北 局 部 山区 
则 偏 小 。 

本 研究 通过 巴基斯坦 气象 局 及 中 国 气象 局 平 
A ,收集 了 印度 河流 域 较为 丰富 的 台 站 实测 资料 。 
经 过 严格 的 气候 界限 值 检 查 和 和 气候 极 值 检 查 后 , 保 
留 20 个 位 于 流域 北部 高 海拔 .复杂 地 形 过 渡 区 的 观 
测 站 (图 1) 。 相 比 再 分 析 资 料 ,本 研究 所 构建 数据 
库 在 印度 河流 域 观测 网 稀少 或 数据 记录 不 完整 的 北 
部 区 域 的 数据 质量 ,得 到 了 一 定 的 提升 。 

基于 格 网 化 的 数据 集 , 本 研究 通过 趋势 分 析 方 
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法 ,揭示 了 整个 流域 气候 变化 的 时 空 差异 性 ,与 前 人 
利用 离散 分 布 的 测 站 资料 得 出 的 气温 降水 时 空 分 布 
规律 基本 一 致 “-”” 。 由 于 流域 一 半 以 上 面积 位 
于 巴基斯坦 境内 ,而且 巴 基 斯 坦 大 部 分 地 区 位 于 印 
度 河流 域 ,已 有 研究 多 数 基于 巴基斯坦 境内 气象 观 
测 资料 分 析 印 度 河 流域 气候 变化 ,忽略 流域 其 他 区 
域 的 气候 变化 特征 。 如 , KHAN 等 中 与 ULLAH 
等 (中 发 现 ,印度 河流 域 巴基斯坦 境内 1960 一 2013 
年 平均 气温 升温 速率 约 为 0.17 ~ 0.33 ©. 
(10 a)”,1980 一 2016 年 降水 量变 率 -10.9 ~12.8 
mm * (10 a) "2*9 ,气温 降水 变 率 均 略 低 于 本 文 结 
果 。 本 研究 发 现 已 有 人 研究 中 被 忽略 的 东北 高 海拔 山 
区 升温 幅度 较 大 ,西部 喀布尔 和 东部 新 德里 降水 增 
幅 较 大 ,这 在 目前 印度 河流 域 气温 降水 研究 结论 中 
未 曾 提 及 。 

印度 河流 域 干旱 年 代 际 波动 与 降水 量 有 着 良好 
的 对 应 关系 。 升 温 背 景 下 ,虽然 年 和 季 降 水 量 在 流 
域 大 部 区 域 都 有 增加 趋势 ,但 是 南部 平原 地 区 ( 
季 ) 潜 在 蒸 散 量 增加 趋势 更 加 显著 ,从 而 导致 流域 
南部 平原 区 表现 出 不 同 程度 的 干旱 加 剧 的 趋势 , 与 
WEN 等 1 的 结论 相同 。 论 文 所 采用 的 标准 化 降水 
Z PATE BX (SPE) 客观 描述 了 历年 印度 河流 域 干旱 ， 
旦 结果 与 标准 化 径流 指数 (SD7) 描述 的 有 关 干 蛙 有 
着 较 好 的 对 应 状态 。 但 个 别 年 份 ,SPE1_6 显示 干旱 
但 SDI 6 却 反映 出 正常 ,说 明 Tronthwaite 方法 获得 
潜在 蒸发 量 对 干旱 略 有 高 估 。 后 续 需 深入 不 同 的 淤 
在 蒸 散 量 理论 在 印度 河流 域 干 旱 检 测 的 比较 研究 。 
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Spatiotemporal distributions of temperature , precipitation ,evapotranspiration , 
and drought in the Indus River Basin 
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2 National Climate Center, Beijing 100081, China; 3 State Key Laboratory of Desert and Oasis Ecology , Xinjiang Institute 
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Abstract: Supported by temperature and precipitation observations from 54 ground-based stations in the Indus 
River Basin and surrounding areas, we used Climate Research Unit temperature and Global Precipitation Climatolo- 
gy Centre precipitation global land surface reanalysis datasets to construct a new gridded monthly dataset covering 
the entire Indus River Basin. This dataset achieved a spatial resolution of 0. 5? x 0. 5? for 1980 —2016. We derived 
potential evapotranspiration values with the Thornthwaite model in order to further analyze the drought characteris- 
tics in the basin according to the Standard Precipitation Evapotranspiration Index. Our results revealed that the areal 
averaged annual mean temperatures increased significantly between 1980 and 2016, with a rate of 0. 30 ^C. - (10 
a) '. The most obvious increase occurred in the early 21st century. On a seasonal timescale, the observed increase 
was much faster in the dry season (November-April; 0.36 °C - (10 a)~') than it was in the wet season ( May- 
October; 0.25 °C + (10 a) ~'). Unlike the monotonic increase observed for temperature, annual precipitation re- 
vealed a decadal oscillation. Driven by the persistent warming, annual and seasonal potential evapotranspiration rates 
have increased significantly. We also detected that drought frequency showed decadal fluctuations similar to those of 
precipitation patterns. Drought was more common in the dry season, but drought intensity is higher in the wet season. 
Meanwhile, the extent of drought-stricken areas experienced a weak increasing trend in the dry season, but a slight 
decrease in the wet season. Spatially, the 1980 -2016 increases in seasonal and annual temperature and potential e- 
vapotranspiration were significant at a 95% confidence level throughout the Indus River Basin, except for a small 
fraction in the northwestern portion of the region. Warming was sharpest in the southern plain and the northeastern 
mountain regions, while potential evapotranspiration rose particularly quickly in the southern plain. Over the same pe- 
riod, annual and wet-season precipitation along the southeast-to-northwest belt from New Delhi to Kabul increased 
significantly , as did dry-season precipitation around Kabul. We detected a significant drying trend during the dry sea- 
son for the southeastern plain and for some parts of the northeastern mountain region. When combined with the drying 
that we saw during the wet season for the northeastern and southwestern mountain regions, our results signal the im- 
portance of strengthening our preparation for climate-related disaster prevention. Effective countermeasures need to be 
identified and implemented in order to mitigate the adverse effects of continued drought intensification. 


Key words: temperature, precipitation; evapotranspiration; drought; spatiotemporal variations; Indus River Basin 


